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In reference[H, Kopnin and Sonin (KS) apply the stan- 
dard BCS model for a two dimensional electron gas with 
the spectrum of Dirac fermions, namely = ap — fi, 
(a ~ ±), where p is the momentum around the Dirac 
point and /i is the chemical potential. Their attempt is 
a generic derivation of superconducting properties dis- 
regarding microscopic details and the sublattice struc- 
ture in graphene. In this comment we argue that apart 
from their derivation of the charge current, the thermo- 
dynamic results of ref. [l[ are not new and were derived 
before in ref.&Q; second, we show that the spectroscopic 
results in ref. |l| such as the superfluid velocity are incon- 
sistent with a Hamiltonian of Dirac fermions. Finally, we 
show that in spite of the fact that their final result for 
the current coincides with the correct result for graphene 
at /i = 0, the derivation for Dirac fermions requires reg- 
ularization, which is only provided by the inclusion of a 
periodic spectrum in the Hamiltonian. 

Ref.[H starts from the usual BCS spectrum for s-wave 
pairing, 



of eigenvectors and also in a different spectrum. 



{^pY + I^P = ±); which is the same 
spectrum derived in Ref. 0, Q from a particular model of 
Dirac fermion superconductivity. The fact that KS do 
not specify a Hamiltonian, however, does not make their 
thermodynamic results more general than those previ- 
ous derivations (as claimed by them) for a trivial reason: 
since the free energy at the mean field level is defined only 
by the spectrum[4], any class of BCS fermionic Hamil- 
tonians which share the same spectrum will have exactly 
the same thermodynamic properties. Since KS start from 
the same BCS spectrum as refjH, they should neces- 
sarily obtain the same results for the gap equation and 
the critical temperature, disregarding any details of the 
matrix structure in the Hamiltonian. KS describe the re- 
sults in Eq. (3) — (11) and the subsequent equation as if 
they corresponded to a new derivation, which is not the 
casejij. 

In the second part of ref. m, KS calculate the super- 
current, j, induced by a uniform flow of the condensate 
with constant momentum ks = Vx, where x — ' r 
is the phase of the superconductor order parameter, 
A = lAje'-^. At the charge neutrality point (/i = 0) 
the Bogoliubov-DcGennes (BdG) equations for a Dirac 
Hamiltonian with s-wave pairing are[6j] 

(p-|-ks)-(7u-f Au = _E{t, -{■p-'ks)-Sv+A*u^ Ev, (1) 

instead of Eq.(2) in ref.[l|, where tr are x,y Pauli 
matrices. These equations result in a different set 



k1 ± 2y/(p • ks)2 + fc2|/\|2^ with distinct spectro- 
scopic properties for finite k^. We note that due to 
particle-hole symmetry, the group velocity of the quasi- 
particles around the Dirac point is zero, whereas the 
particle- hole charge current is finite [7]. This symme- 
try argument shows that the spectrum derived in ref.jlj 
(which gives a finite superfluid velocity at half filling) 
is inconsistent with any BdG Hamiltonian of Dirac 
fermions, and therefore is not applicable to graphene. 

Finally, using a covariant momentum in the BdG 
Hamiltonian, H , namely k^ = Vx — A, (or ks = —A 
by a gauge choice), where A is the vector potential, 
the current follows from j — —d{H)/dA. From Eq. 
(1) one finds j oc {D - 2| A|tanh[|A|/(2T)]}A, where 
D ^ |A| is the band width, which accounts for the or- 
bital paramagnetic response of the lower band electrons, 
overwhelming the diamagnetism. This term is absent 
from the current definition of[lj without justification[§]. 
The diamagnetism is recovered only if one includes the 
full spectrum, = a|0p| = o.\(t>'p + *0pli where 0p 
is a periodic function. In that case, the graphene 
BdG equations, i(l>p-A^x + 4>'-L- /<<^v)u + Aw = Eii, and 

-{4>'p+A(^x+4>'l,+j^(Ty)V' 



"p-K'-'Vl 

A*u = Ev, will give ji oc {Si 



|A|2Ep|5p.0pl'i^p'5Bjtanh(£;p/2T)/i?p]}A„ (i =x,y 
directions) where Si is a surface term which is regularized 
by the Brillouin zoneQ. 
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2 



metal) case, where j = 0. In contrast with metals, the 
free energy for Dirac fermions depends on _D [^j, and the 
Ginzburg-Landau current derived from it requires regular- 
ization. 



[9] Si = —^piydp-{a{dp-Ep)n{aEp)} is exactly zero when 
calculated in the whole Brillouin zone, and is oc D for a 
linear (non-periodic) dispersion. 



